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1 Introduction 



Theories of extra spatial dimensions, in which the observed universe is reahzed as a brane 
embedded in a higher dimensional bulk, have attracted a lot of attention in the last few 
years. In this viewpoint, ordinary matters and gauge fields are trapped on the brane but 
gravitation and possibly non-standard matter can propagate through the entire spacetime 
[1,2,3]. In a cosmological perspective, braneworld scenarios have the capability to explain 
some of the new achievements of observational cosmology such as late-time positively 
accelerated expansion. Although some of these models predict deviations from the usual 
4-dimensional gravity at short distances, the model proposed by Dvali, Gabadadze and 
Porrati (DGP) [1] is different in this regard since it predicts deviations from the standard 
4-dimensional gravity over large distances. In fact, DGP scenario is infra red modification 
of general relativity. 

On the other hand, impacts of Lorentz invariance violation (LIV) on cosmology have 
been studied by some authors [4,5]. For instance, this issue has been studied in the 
context of scalar- vector-tensor theories [4]. It has been shown that Lorentz violating 
vector fields affect the dynamics of the inflationary models. One of the interesting feature 
of this scenario is the fact that exact Lorentz violating inflationary solutions are related to 
the absence of the inflaton potential. In this case, the inflation is completely associated 
with the Lorentz violation and depends on the value of the coupling parameters [5]. 
One important observation has been made recently in references [6] which accelerated 
expansion and crossing of the phantom divide line with one minimally coupled scalar field 
in the presence of a Lorentz invariance violating vector field has been shown as a result 
of interactive nature of the model. We know, an important consequence in quintessence 
model is the fact that a single minimally coupled scalar field is not suitable to explain 
crossing of the phantom divide line, a; = —1 [7]. However, a single but non- minimally 
coupled scalar field is enough to cross the phantom divide line by its equation of state 
parameter [8] . It has been shown that within an interactive picture, a minimally coupled 
scalar field in the presence of a Lorentz violating vector field can evolve to take phantom 
phase [6]. This happens due to wider parameter space prepared by inclusion of Lorentz 
invariance violating vector field and its interaction with scalar field. 

From another viewpoint, currently it is well-known ( based on various observational 
data) that our universe has entered the stage of a positively accelerate expansion around 
the redshift z < 1; see [9] and references therein. The equation of state (EoS) parameter u 
responsible for the acceleration of the Universe has been constrained to be close to = — 1. 
In this regard, the analysis of the properties of dark energy from recent observational 
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data mildly favor models of dark energy with uj crossing —1 line in the near past. So, the 
phantom phase EoS with a; < — 1 is still mildly allowed by observations. Recently, there 
have been a number of attempts to realize the phantom phase EoS, the simplest model 
which realizes the phantom EoS is provided by a nonminimally coupled scalar field [7]. 
Other examples are in the spirit of braneworld models of dark energy. In these models, 
crossing of the phantom divide line and late-time acceleration are studied extensively ( 
see for instance [10,11,12]). 

With these preliminaries, construction of yet another theoretical framework which 
combines braneworld effects and Lorentz invariance violation to realize positively accel- 
erate expansion and transition to phantom phase is an interesting challenge. With this 
motivation, in which follows, we construct a new dark energy model to realize crossing of 
the phantom divide line and explanation of other observational achievements such as late 
time acceleration. In this regard, we study cosmological dynamics of a Lorentz violating 
DGP-inspired braneworld scenario. By implementing local Lorentz violation in a grav- 
itational setting due to the existence of a tensor field with a non-vanishing expectation 
value, and then coupling of this tensor field to gravitational sector and matter (a scalar 
field) , we study late-time acceleration and transition to phantom phase of the scalar field. 
The simplest example of this approach is to consider a single timelike vector field with 
fixed norm. This vector field picks out a preferred frame at each point in space-time and 
any matter field coupled to it will experience a violation of local Lorentz invariance. A 
special case of this theory was firstly introduced as a mechanism for Lorentz- violation by 
Kostelecky and Samuel in Ref. [13]. In curved spacetime, however, there is no natural 
generahzation of the notion of a constant vector field (since V ^u^ — generically has no 
solutions); we must therefore allow the vector field to have dynamics, and fix its norm 
by choosing an appropriate action for the field. We will show that as a result of Lorentz 
invariance violation, there is an interacting term in the dynamics of scalar field which 
affects cosmological dynamics of the model considerably and is responsible for transition 
to phantom phase of the scalar field. 

Nevertheless, one point should be stressed here: we know that DGP braneworld sce- 
nario explains accelerated expansion of the universe via leakage of gravity to extra di- 
mension [14]. But in this scenario the EoS parameter of dark energy never crosses the 
oj{z) — —1 line and universe eventually turns out to be de Sitter phase. On the other hand, 
in this setup by incorporating a single scalar field (ordinary or phantom) on the brane, 
one can show that EoS parameter of dark energy crosses the phantom divide line [15]. 
The question then arises: why we need to further generalization of this braneworld setup? 
One important reason lies in the fact that DGP setup suffers from ghost instabilities and 
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it is important to add new ingredients to original setup to overcome this shortcoming. 
We thing Lorentz invariance violation may help us to construct a wider parameter space 
with potential to overcome this problem. 

We begin by an overview of the basic equations of motion for the most general theory 
of a fixed-norm vector field with a generalized action S — SBuik+S Brane where S Brane = 
Seh + Stf, + Sjn + Su- Then we generalize this setting to the model universe with DGP- 
inspired action and we investigate cosmological consequences of this setup. 



2 A Lorentz Violating DGP-Inspired Braneworld Sce- 
nario 

As it is well known, brancworlds arc often studied within the framework of the 5D Einstein 
field equations projected onto the 4D brane [16]. To study impact of Lorentz invariance 
violation on the cosmological dynamics of DGP setup, we consider a vector field along 
with the extra dimension. So, a local frame at a point in space-time is inevitably selected 
as the preferred frame. In other words, the existence of the brane defines a preferred 
direction in the bulk. We study the effects of local Lorentz violation on the dynamics of 
the brane by inclusion of this vector field in the action. This additional field modifies the 
AD Einstein equations with cosmological implications which we investigate by studying 
the resulting Priedmann equation on the brane. 

The action of the Lorentz violating DGP scenario in the presence of a minimally 
coupled scaler field and a vector field on the brane can be written as the sum of two 
distinct parts 

S — SbuIH + ^Branej (l) 

where Ssuik and Ssrane = Seh + Ss + + Sy, are defined as follows 
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Ssuik = / d'x'^^n, (2) 



s 



Brane 



/ rji^ 1 



Jj/=0 

(3) 

where y is coordinate of the fifth dimension and we assume brane is located at y = 0. m\ 
and m\ are fundamental scales in the bulk and brane respectively, qab is five dimensional 
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bulk metric with Ricci scalar 7?., while q^i, is induced metric on the brane with induced 
Ricci scalar R. qab and q^j^t, are related via q^j^t, — 5i/^5i,^ qab- K is trace of the mean 
extrinsic curvature of the brane defined as 



K 



— lim ( 
2 6^0 V 



+ 



y=+e 



(4) 



and corresponding term in the action is York-Gibbons-Hawking term [17]. This action is 
allowed to contain any non-gravitational degrees of freedom in the framework of Lorentz 
violating scalar- vector-tensor theory of gravity. As usual, we assume uf^Ufj, — —1 and that 
the expectation value of vector field is < 0\u^Uij,\0 >— —1 [18]. {i — 1, 2, 3, 4) are 

arbitrary parameters with dimension of mass squared and A is a Lagrange multiplier. Note 
that are mass scale of Lorentz symmetry breakdown [4,18,19]. In which follows, we 
neglect quartic self-interaction term, u'^u'^V ^u'^V ^Ua- Some cosmological consequences 
of this term in the action are studied in Refs. [4,19] for ordinary AD framework. In this 
setup, the preferred frame is selected through the constrained vector field u'^ and this 
leads to the violation of the Lorentz symmetry. 

The ordinary matter part of the action is shown by Lagrangian = ^miQuu,^^) where 
ip is matter field and corresponding energy-momentum tensor is 

5C 



T, 



-2 



(5) 



The pure scalar field Lagrangian, 
energy-momentum tensor 



Sq^"" 

= -i^M'^V^^V^t/) - V{(f)), yields the following 



T, 



l^(lAV(t>f-q^r.V{4>). 



(6) 



The energy-momentum tensor of the Lorentz violating vector field is defined as usual by 



rp{u) 



5q' 



(u) 



(7) 



The Bulk- brane Einstein's equations calculated from action (1) are given by 

ml (Uab - \gABn) + m^^^^^B^ [r^, - 5{y) = Sa'^Sb'^T ^J{y) (8) 

where T^,^ = T^,^ -|- r^^ -|- Tj^ . From equation (8) we find 



G 



and 



AB 



Gfj,u — 



T^AB — -^9abT^ 



R 



T 



(9) 



(10) 
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as Einstein's equations in the bulk and brane respectively. The corresponding junction 
conditions relating the extrinsic curvature of the brane to its energy-momentum tensor, 
have the following form ( see [20] and related references therein) 



lim 



y=-e 



1 



T 



y=o m% 



y=0 



(11) 



We start with the following line element to derive cosmological implications of our model 

ds'^ = q^^dx'^dx'' + t)dy'^ = -n^iv, t)dt^ + a^{y, t)-fijdx'dx^ + b^{y, t)dy^. (12) 

For such a metric, to solve Einstein equations in the presence of a fixed-norm vector 
field, the vector field must respect spatial isotropy, at least in the background (though 
perturbations will generically break the symmetry). Thus the only component that the 
vector can possess is the timelike component. Therefore, we take the constraint u'^ ~ 
0, 0, 0) where jV is a lapse function. After performing required algebra, we set n(0, t) = 
1 and jV' = 1 in which follows. The scale of the universe is determined by a{y, t) and 
is a maximally symmetric 3-dimensional metric defined as 



X j^X j 



(13) 



where k = —1,0, 1 parameterizes the spatial curvature and = XiX^. Wc assume that 
scalar field depends only on the proper cosmic time of the brane and we adopt the gauge 
h'^{y,t) = 1 in Gaussian normal coordinates. 
Now total energy density and pressure are given as follows 



Ptot ^ Pm + Pu + Pel, 



(14) 



and 



Ptot = Pm + Pu + P4>- (15) 

Here we assume that energy-momentum tensor of ordinary matter on the brane has a 
perfect fiuid form with energy density pm and pressure pm so that T^j, = {pm+Pm)Nfj_N^ + 
PrnQiiu where Nfj_ is a unit timelike vector field representing the fiuid four- velocity We also 
assume a linear isothermal equation of state for the fiuid Pm = (7m~l)Pm that 1 < 7m ^ 2. 
Energy density and pressure of minimally coupled scalar field are given as follows 

1 



and 



Pep 



P<t> 



.2v? 



+ n'V{(t>) 



02 _ V{ct>) 



y=o 



(16) 



(17) 
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where a dot denotes the derivative with respect to cosmic time t. The stress-energy for 
the vector field also takes the form of a perfect fiuid, with energy density given by [4,5] 



Pu 



(18) 



and pressure 



p„ = /3//2[3 + 2-^ + 2 



/3 



HP 



(19) 



where we have defined the parameter p as P = {Pi + 'iP2 + P3) and H = is the Hubble 
parameter. In the absence of vector field, that is, when all Pi = 0, the above equations 
reduce to the conventional ones and in the case P = const., the above equations arc lead 
to the equations given in [19]. For future references, the total equation of state parameter 



defined as ujtot = — or 



^tot 



Vm+Pu+ . 
Pm + Pu + 



takes the following form 



rp^ + PH' [3 + 2^ + 2^1 + - V{<P) 



^tot 



HJ3_ 
1^2 



y=o 



Pm-3PH^ + icj>^ + n^V{cf>) 
Now, the effective Einstein equations on the brane are given as follows [21] 

n, 



(20) 



Gfxu — 



ml 



where 



and 



n 



= CijKL e>' e^'g^ 9 



L 



[21] 



[22] 



(23) 



where Cjjkl is five dimensional Weyl tensor and is the spacelike unit vector normal 
to the brane. Using equation (21) we find 



n% 

ml 



cQ 
^ 



where for FRW universe we have 



Similarly, for space components we have 



ml 



(24) 



(25) 



(26) 
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where 

G'j = -(2ij + 3H^ + ^)(5V (27) 

Now, using equation (22) we find n% = and = - j^T% (T% - 2Tii)5*j. 

Also, the time and space components of the total energy- momentum tensor are given by 

T\^-pt,t-mlG\ (28) 

and 

T^^-ptot^i-m^G^,-, (29) 

where ptot and ptot are given by (14) and (15). These equations lead us to the following 
effective Priedmann equation on the brane 



(30) 



This equation contains Lorentz invariance violation via existence of j3. One can easily 
show the possibility of existence of two different branches of this DGP-inspired Fried- 
mann equation. On the other hand, as the tensor structure in the action is very different 
from the original DGP model, we can expect the above action contains no ghost. In fact, 
we have many parameters here, hence there is a chance to remove the ghost instabilities 
from our model. As we will show, by suitable fine-tuning and constraining parameters 
of this model with observational data, this model has the capability to explain latc-timc 
acceleration and other cosmologically interesting issues. For future reference, we note that 
using Codazzi equation one can show that V^S^i, = [21]. 



3 Cosmological Aspects of the Model 

Late-time positively accelerated expansion and transition to phantom phase of scalar 
field as a dark energy component are two interesting challenges of modern cosmology. 
Therefore, in which follows, we study late-time acceleration and dynamics of equation 
of state parameter, a;<^(i) in this Lorentz violating DGP setup. The equation of state 
parameter of scalar field on the brane is given by 



P4' 



102 + ri^v{(l)) 



(31) 

y=o ^ ' 
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where from now on we set n(0, t) — 1. Dynamics of u!^{t) can be obtained in two different 
viewpoints. Firstly, Priedmann equation (30) in the absence of ordinary matter on the 
brane and with flat spatial geometry (A; = 0), can be rewritten as follows 



1 • l2 



We take E'^o + 'iHE'^o — which an integration gives S^o = ^ where Sq is integration 
constant [21]. Dynamics of scalar field can be deduced from this equation directly 



0' = 6H\P + ml) - 2V{(f)) + 12mle\ H^ + ^ (33) 

y o(X 

where e = ±1, corresponding to two possible embedding of DGP braneworld. To deter- 
mine exact dynamics of scalar field, we need to specify functional form of potential, V{(f)). 
In this framework, we use two well-known potentials: V{(j)) = and an exponential 
potential as V{<j)) = ^exp - ^^0) [22]. 

Secondly, we see that equation (33) depends on the potential of scalar field. With this 
facility, we can obtain a potential-independent equation of dynamics for scalar field. The 
energy equation for vector field u is 

Pu + 3H{p^+Pu)^+3H''p (34) 

and for the scalar field we find 

+ 3H{p^ + p^) = -3H^p. (35) 

There is a non-conservation scheme in this setup due to energy-momentum transfer be- 
tween scalar and vector fields. This is very similar to the case studied by Zimdahl et al 
[23] . As they have shown, a coupling between a quintessence scalar field and a cold dark 
matter (CDM) fiuid leads to a stable, constant ratio for the energy densities of both com- 
ponent compatible with a power law accelerated cosmic expansion. In fact this coupling 
is responsible for accelerated expansion and possible crossing of phantom divide line. In 
our Lorentz invariance violating scenario this coupling is present between scalar field and 
vector field leading to an interactive picture. Note that one can consider also an interac- 
tion between ordinary matter and vector field on the brane. Very recently, authors of Ref. 
[24] have investigated the cosmological evolution of an interacting scalar field model in 
which the scalar field has an interaction with the background matter via Lorentz violation 
in 4-dimensional model. 

Nevertheless, the total energy in the presence of both scalar and vector fields is conserved 

p + 3H(p + p)^0, (p^p^ + p^). (36) 
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This energy conservation equation can also be obtained by equating the covariant diver- 
gence of the total energy-momentum tensor to zero, since the covariant divergence of the 
Einstein tensor is zero by its geometric construction. It follows from contraction of the 
geometric Bianchi identity. 

We obtain dynamics of the scalar field by differentiating equation (16) with respect to t 
and then using equation (35) to find 

4> + 3H^ + 3H^P,^ + V^ = 0. (37) 

By differentiating equation (32) with respect to t and using equation (37) we have 

= T2m^(^ + -i^)^ - - 2(5H, - (38) 

where we assume H and j3 arc depended on (j) in forthcoming arguments and for simplicity, 
wc set £q = 0. By substituting equation (38) into the Friedmann equation (32), the 
potential of the scalar field in our model takes one of the following forms 

, ^, 1 12H + 6m3'//'m4' - + 6HVP 
^^'^^ = 2 

and 

1 -12// + Gma'// W - 0W + 
^^'^^ = 2 

where can be obtained from equation (38). We emphasize that equation (38) has several 
solutions, here we just consider one real root of this equation to proceed further. On the 
other hand, equation governing on dynamics of itself has no analytical solution. So, one 
can try to find some intuition by numerical analysis of the parameter space of the model. 

To obtain dynamics of equation of state parameter, we note that there are two gov- 
erning equations on the dynamics of //^ as follows 



^ ^3 2pm3^ + + m3^ ^^^^ 

and 

2 1 P^ma^ + PPci> + 6m/ - 2 yJSmi^p^ms^ + 3m4^P p^ + 9m4^^ 
^ ^3 2pms^ + p^ + ms^ ^^^^ 

Using these two branches of the model and also equation (35), we find 

Ei + E2(l+u;^)^-PEs (43) 
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where we have defined 



H (3 HrriiC' + 13 H m/" + m^'^H niA 



H 



+ 



(44) 



3 ei 777,4^ + €2 



E. 



■014'' +11 - 1113-)' nij '' +2 Ucj ii:i-riii['r) 
3 



7714'^ 



and 



2 €3 ^77743 (77749 + H^(P + m^^f + 2 He4 {(3 + mz'^))mi^ 

H (7774=^ {(3 + 7773") + 2 £4) (/3 + 7773^)' 

(/5 + ^3^)^ 7774^ + 2 He^ {(3 + 7773^) + 2 7774^ 



(45) 



(46) 



(77743 {(3 + m.J) iJ + 2 €4) (/5 + 77732)' 
In these relations we have defined ei = £2 = £3 = €4 = ±1 where have been appeared due 
to algebraic structure of the model. Although all of these quantities have the same value, 
we have to save them in forthcoming equations because of several permutations of signs in 
equations of cosmological dynamics. In fact, as we will show, suitable and simultaneous 
choices of these quantities are important in analysis of parameters space. Now, equation 
of state parameter takes the following complicate form 



1 



- 6 62 ((3 + 7773') Hei (1/2 7774^ (/3 + 7773') + 64) X 



\ 



7774^ + H^(3 + 7773^)^ 7774^ + 2 ^£4 {(3 + 7773^) 
7774^ 



2 P,^ <P 63 ^^7743 (77749 + H^p + 7773^)' 77743 + 2 ii'£4 iP + 77732))7774'- 

37774^ (/? + 7773')' (2/3 + £1 7774^) H''- 
6 ( (1/3 H,^ 4> 7773' + 2/3 H^^ <j)Pms'' + e, £4 + 1/3 P^H^ <j)) m/+l/3 P,^ £4) {p + 7773') H- 



2 (7774 V,,/. + £4 [P + 7773')') 



7774 + £2 



77749 + H^p + 7773^)2 77743 + 2 i/£4 {P + 77732) 



7774" 



X (/? + 7773') ' H-hi-^ (1/2 7774' {p + 7773') ii" + £4) 



(47) 
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This equation will be used to perform numerical analysis of the model. As another 
important cosmological parameter, the deceleration parameter q which is defined as 




(48) 



using equation (32), takes the following form 



(-12/3^ 



- -2^665^ (6m^m^ + {2Hp + mg^) (^-3/3H^ + V(0) + 00) )H^m4^+ 

12m4^) H^-12P (ms^ - l/2p) i/3+(-3m3^ + 3 pm^^) H'^-2p (v{(j)) + 00) H 
-ma' (V^(0) + 00) [(2H(3 + m3^y H^]~\ (49) 



Recent observations of distant type la supernovae and other observational data [9] indi- 
cate that q is currently negative; that means the expansion of the universe is positively 
accelerated. This is an indication that the gravitational attraction of matter, on the cos- 
mological scale, is more than contraction by negative pressure dark energy in the form of 
quintessence. In which follows, we study cosmological implications of our model focusing 
on late-time acceleration and scalar field dynamics as a candidate of dark energy and pos- 
sible transition to phantom phase. We show that the Lorentz symmetry breaking scalar 
field can be treated as a good candidate for the role of the dark energy source. 
Before proceeding further, we should address some important and related issues here. 
Firstly, one should be careful to choose the appropriate equation of state for components 
that are used to describe the universe energy- momentum content. As we have emphasized 
earlier, a suitable coupling between a quintessence scalar field and other matter content 
can leads to a constant ratio of the energy densities of both components which are com- 
patible with an accelerated expansion of the universe or crossing of the phantom divide 
fine (for more details see [23] and reference therein). In this respect and for instance, the 
holographic dark energy models studied in Ref. [7] have the phantom phase by adopt- 
ing a native equation of state, whereas the authors in [25] have found accelerating phase 
only using the effective equation of state. Based on these arguments, we should explain 
what kind of equation of state is used for observing the nature of mixed fluids here. In 
our model, we have three sources of energy-momentum: 1- standard ordinary matter, 2- 
scalar fleld as a candidate of dark energy and 3- energy-momentum content depended 
on Lorentz violating vector fleld. Here we assume that standard matter has negligible 
contribution on the total energy-momentum content of the universe and we can consider 
a constant linear isothermal equation of state as Pm — {'Jm — l)Pm that 1 < 7m < 2 for 
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it. For other two energy- momentum contents, it is possible to use the "trigger mecha- 
nism" to explain dynamical equation of state [26]. This means that we assume scalar- 
vector-tensor theory containing Lorentz invariance violation which acts like the hybrid 
inflation models. In this situation, vector and scaler field play the roles of inflaton and 
the "waterfall" field respectively [26]. In this regard, we can fine-tune parameter m and 
other parameters to obtain best fit model using the observational data. Of course, an 
attractor solutions and fine-tuning in Lorentz violation model for suitable inflation phase 
has been studied in Ref. [5]. Therefore it is reasonable to expect that one of them will 
eventually dominate to explain inflation or accelerating phase and crossing of phantom 
divide line. We should emphasize that the model studied in this paper belongs to a wider 
class of Lorentz- violating theories exhibiting the phantom behavior (see for instance Refs. 
[4,5,27] for a number of Lorentz- violating models). One can extend this framework to 
understand issues such as transient phantom stage, super-horizon ghosts and to deal with 
the question that how generic are the features obtained in this particular model of late- 
time de Sitter attractor. Some of these issues have been discussed in Ref. [28]. One more 
direction in this framework is to modify our model in such a way that it find the capability 
to describe inflationary epoch rather than the late-time acceleration. As has been pointed 
in [28] , this type of model may give rise to some distinct features in the CMB spectrum. 



4 Numerical Analysis of the Parameters Space 

In this section we study late-time acceleration and possible transition to phantom phase 
of the scalar field. Incorporation of the Lorentz invariance violation in the model provides 
a wider parameter space (relative to the case with one Lorentz preserving scalar field) 
which leads to more suitable framework for explanation of these interesting cosmological 
aspects. To show the validity of this statement, we need to solve equations (47) and (49). 
In the first stage, we obtain dynamics of scalar field using equation (33). This goal will 
be achieved only if the Hubble parameter H{(j){t)) and the vector field coupfing, /3{(p{t)) 
are known a priori. In which follows, our strategy is to choose some suitable and natural 
candidates for the Hubble parameter H{(p{t)) and the vector field coupfing P{(p{t)). Then 
we focus on possible crossing of the phantom divide line and realization of the universe 
late-time acceleration. We obtain suitable domains of parameters space which admit late- 
time acceleration and crossing of the phantom divide line by equation of state parameter. 
Probably this Lorentz violating DGP-inspired model has some important consequences in 
the spirit of cosmology ( such as possible realization of bouncing solutions) and particle 
physics, but here we focus only on late-time acceleration and transition to phantom phase 
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of the scalar field. 

We consider a general case where both the vector field couphng and the Hubble pa- 
rameter are functions of scalar field defined as follows 

H = Ha(t)^ , = m0« (50) 

where Hq and m are positive and constant parameters. Note that scalar field itself is 
depended on the cosmic time, t. In this case, dynamics of scalar field with potential of 
the type V{(j)) — X'(f)^ is given by 

<j) = 6//o^0^^(m(/)^ + ms^) - AV^ + 12m4^ei/o0^ (51) 

and for potential of the type V{4>) = Vocxp ^ — ^J^^'pj , we find 



= 6//o'0'^ {m(f>^ + ma') - exp f - ^ + 12m/ei/o</'^- (52) 

The question then arises here: what choices of space parameters lead to analytical so- 
lutions of the equations (51) and (52)? The answer to this question is summarized in 
tables 1 and 2. We note that as table 2 shows, for exponential potential there is very 
limited possibility to obtain analytical solution for scalar field dynamics. Nevertheless, 
we can find analytical solution in some especial cases. For instance, form equation (51) 
with ^ = — 1 and C = 1 we find the following analytical solution for dynamics of scalar 
filed 



-{t+Ao)^JQHo^m.3^^-\' 



(6i/oW-A')'^']"' (53) 



where is an integration constant. We emphasize that this analytical solution of equa- 
tion (51) is obtained under some especial choices of parameter space. Without these 
choices it is impossible to find closed analytical solution for 0(i). We use this solution 
for our forthcoming numerical analysis. Figure 1 shows dynamics of as described by 
equation (53) for two different cases. 

Now by substituting equation (53) into equation (49), we can analyze dynamics of 
deceleration parameter q with respect to cosmic time t. This has been shown in figures 2 
and 3 and corresponding results are summarized in table 3. 
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Figure 1: Dynamics of scalar field in two especial cases. 



Table 1: Acceptable range of ^ and ( to have analytical solution of the scalar field equation 

(51) 



values of ( and ( 


aualytical solutiou for e(iuati()u (51)? 


^ ^ C) aud ( < C) 


uo 


{ = aud ( = C) 


yes 


^ = 1 and C = 


yes 


^ = 2 and C = 


yes 


^ > 1 and C > 1 


no 


e = -1 and C = 1 


yes 


^ = -2 and C = 1 


yes 


other ^ and ( 


should be examined 
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Table 2: Acceptable range of ^ and C to have analytical solution of the scalar field equation 
(52) 



values of C, and ^ 


analytical solution for equation (52)? 


e 7^ and C ^ 


no 


{ = and C = 


yes 


C = 1 and C = 


no 




Figure 2: Variation of deceleration parameter q relative to cosmic time t for e = +1 and = ±1 
with scalar field potential of the type V{(f)) = X'cfP'. 



Table 3: Summary of results from figures 2 and 3. 



e 


£5 


deceleration parameter q 


late-time acceleration ? 


+1 


+ 1 


negative 


yes 


+ 1 


-1 


negative 


yes 


-1 


+1 


negative 


yes 


-1 


-1 


negative but almost constant 


yes 
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Figure 3: Variation of deceleration parameter q relative to cosmic time t for e = — 1 and £5 = ±1 
with scalar field potential of the type V^(^) = X'cfP'. 



In which follows, we consider equation of state as given by equation (47) with equation 
(53) for dynamics of scalar filed in order to study crossing of the phantom divide line in 
this setup. For this purpose, we set ^ = —1 and C = 1 to be more specific. The results 
of numerical calculations are shown in figures 4, 5, 6 and 7. All of these figures show 
the possibility of transition to the phantom phase of the scalar field. It is interesting to 
note that as these figures show, transition from quintessence to phantom phase and from 
phantom phase to quintessence is possible in this braneworld scenario with appropriate 
choices of model parameters. Since one minimally coupled scalar field in the presence of 
Lorentz invariance symmetry cannot realize phantom divide line crossing, observation of 
this crossing in our model is a result of interactive nature of the model due to Lorentz 
invariance violation. This is manifested in our model via existence of the term such as 
— 3/9if^ in Friedmann equation (32) or —3$H'^ in equation (35). 

We should stress here that by neglecting ordinary matter content { Pm — — Pm) in 
equation (20), we find 



^tot 



(54) 



The same procedure as described above, leads us to the dynamics of cUfot ^ shown in 
figure 8. We see that for e = +1, cotot never crosses the phantom divide hne. 
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Figure 4: Crossing of phantom divide line. 
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Figure 5: Crossing the phantom divide line. 
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Figure 6: Crossing the phantom divide line. 
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Figure 7: Crossing the phantom divide line. 
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Figure 8: Dynamics of LOtot- There is no crossing for e = +1 branch of the model. 

5 Summary and Discussion 

In this paper we have studied a DGP-inspired braneworld scenario where the idea of 
Lorentz invariance violation has been incorporated by specifying a preferred frame through 
the introduction of a dynamical vector field normal to our brane. The Einstein field equa- 
tions obtained on the brane are modified by additional terms emanating from the presence 
of the vector field. This model breaks the AD Lorentz invariance in the gravitational sec- 
tor. As a new mechanism for crossing of the phantom divide line by equation of state 
parameter, we have shown that by a suitable choice of parameters of the model, it is 
possible to have phantom divide line crossing and realizing late-time acceleration in this 
Lorentz invariance violating context. 

More importantly, we need to show the stability of the solutions in this Lorentz vi- 
olating DGP-inspired model. It has been shown that the self-accelerating branch of the 
DGP model contains a ghost at the linearized level [29]. The ghost carries negative en- 
ergy density and it leads to the instability of the spacetime. The presence of the ghost 
can be related to the infinite volume of the extra dimension in the DGP setup. When 
there are ghost instabilities in the self-accelerating branch, it is natural to ask what are 
the results of solution decay. One possible answer to this question is as follows: since 
the normal branch solutions are ghost-free, one could think that the self-accelerating so- 
lutions may decay into the normal branch solutions. In fact for a given brane tension, 
the Hubble parameter in the sclf-accclcrating universe is larger than that of the normal 
branch solutions. Then it is possible to have nucleation of bubbles of the normal branch 
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in the environment of the self- accelerating branch solution. This is similar to the false 
vacuum decay in de Sitter spacetime. However, there are arguments against this kind of 
reasoning which suggest that the self- accelerating branch does not decay into the normal 
branch by forming normal branch bubbles [29]. It was also shown that the introduction 
of a GaussBonnet term for the bulk does not help one to overcome this problem [30] . In 
fact, it is still unclear what the end state of the ghost instability is in the self-accelerated 
branch of DGP inspired setups (for more details see [29]). On the other hand, it seems 
that introduction of Lorentz violating vector field on the branc provides a new degree 
of freedom which may provide a suitable basis to treat ghost instability. Field theories 
of vector fields are very restrictive and gauge invariant Lagrangians are the ones that 
guarantee that the zeroth component of the gauge field is not propagating and, therefore 
there are no ghosts. However, in the action of our model as given by equation (3), the 
brane Lagrangian for the vector field is not gauge invariant. Therefore, it is possible that 
this vector field action by itself can introduce a ghost, instead of curing the one of gravity. 
This is not necessarily the case in the presence of vector field condensates [31,32]. On the 
other hand, for certain timelike vector theories with spontaneous Lorentz violation, suit- 
able restrictions of the initial-value solutions are identified that yield ghost-free models 
with a positive Hamiltonian[33]. This restrictions can be imposed on parameters such as 
Pi in our model ( see also [19,30]). Therefore, in our Lorentz invariance violating setup 
there is the capability to overcome instabilities both in gravitational and vector field sec- 
tors. Anyway, if the Lorentz violation annihilate the ghost, this will be a great progress in 
the cosmology. We think that our Lorentz invariance violating model on the DGP brane 
has the capability to solve ghosts instabilities problem due to its wider parameter space 
relative to other existing scenarios. 

Within a similar viewpoint, very recently Zen et al [24] have studied the cosmological 
evolution of an interacting scalar field model in which the scalar field has its interaction 
with dark matter, radiation, and baryon via Lorentz violation in 4£) standard model. They 
proposed a model of interaction through the effective couphng parameter, /9, — — 
They also determined all critical points and studied their stability of Lorentz violation 
model. On the other hand, Mariz et al [34] studied the fact that the Lorentz-breaking 
parameter can be treated as a natural explanation of the extremely small value of the 
cosmological constant. Thus, the Lorentz symmetry breaking introduces a mechanism 
for the arisal of a non-zero but very small cosmological constant, and therefore providing 
an acceptable solution for cosmological constant problem. We believe that Lorentz sym- 
metry breaking fields can be treated as a good candidate for the role of the dark energy 
source and these types of models have the capability to address other cosmological issues 
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with better adaptability than other models. Finally we should stress on some open issues 
in this field: one issue concerns that we have provided a scenario of a scalar field with 
Lorentz violation to realize the crossing of the phantom divide. One maybe curious to 
know whether this scenario can make some interesting predictions for observations. For 
example the interactive terms in this model is very important, but can this interaction 
affect the formation of large scale structure, or the occupation of dark matter? This is 
a very interesting issue, since as we know large scale structure and the measure of dark 
matter are usually strongly constrained by CMB observations. Moreover, we think this 
model is able to give a bouncing solution of the universe. The reason lies in the fact that 
as has been shown in Ref . [35] , any models with crossing of phantom divide line have the 
capability to realize a bouncing solution. 
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